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Abstract 

We report new constraints on extra-dimensional models and other physics beyond the Standard 
Model based on measurements of the Casimir force between two dissimilar metals for separations 
in the range 0.2-1.2 //m. The Casimir force between an Au-coated sphere and a Cu-coated plate 
of a microelectromechanical torsional oscillator was measured statically with an absolute error of 
0.3 pN. In addition, the Casimir pressure between two parallel plates was determined dynamically 
with an absolute error of 0.6 mPa. Within the limits of experimental and theoretical errors, 
the results are in agreement with a theory that takes into account the finite conductivity and 
roughness of the two metals. The level of agreement between experiment and theory was then used 
to set limits on the predictions of extra-dimensional physics and thermal quantum field theory. It 
is shown that two theoretical approaches to the thermal Casimir force which predict effects linear 
in temperture are ruled out by these experiments. Finally, constraints on Yukawa corrections to 
Newton's law of gravity are strengthened by more than an order of magnitude in the range 56 nm 
to 330 nm. 

PACS numbers: 03.70.+k, 12.20.Fv, 12.20.Ds, 42.50.Lc 
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I. INTRODUCTION 



Many extensions of the Standard Model, including supergravity and string theory, exploit 
the Kaluza-Klein idea that the true dimensionality of space-time is N = 4 + n, where the 
additional n spatial dimensions are compactified at some small length scale. For a long time 
it was generally believed that the compactification scale was on the order of the Planck 
length I pi = \[G ~ 10 _33 cm, where G is the Newtonian gravitational constant, and units 
are chosen so that h = c = 1. The corresponding energy scale Mpi = 1/y/G ~ 10 19 GeV is 
so high that direct experimental observation of the effects of extra dimensions would seem 
impossible at any time in the foreseeable future. 

The situation changed dramatically with the proposal of models for which the com- 
pactification energy may be as low as the extra-dimensional Planck energy scale, Mpf = 
l/G l ^ 2+n \ which is assumed to be of order ITeV 

HQ. 

(Here is the fundamental 

gravitational constant in the extended iV-dimensional space-time). Note that this proposal 
eliminates the hierarchy problem since the characteristic energy scales of gravitational and 
gauge interactions coincide. In order to be consistent with observations, the usual gauge 
fields of the Standard Model are presumed to exist on 4-dimensional branes whereas gravity 
alone propagates into the iV-dimensional bulk. 

Constraints on these new lower energy scale compactification models can be obtained by 
investigating their predictions in accelerator experiments 0,0,0], astrophysics 0, Q, 0], and 
cos.no.ogy 0130. More modeMndependen. .imits, however, can be ot.ta.ned from tests 
of Newtonian gravity. In extra-dimensional models with large compact extra dimensions , 
the Newtonian gravitational potential acting between two point masses acq uires a Yukawa 
correction for separations much larger than the compactification scale while for 

models with non-compact but warped extra dimensions, the corrections are power-laws 2]. 
For two interacting macroscopic bodies, either of these corrections would give rise to a new 
(so called "fifth") force coexisting with the usual Newtonian gravitational force and other 
conventional Standard Model interactions, such as Casimir and van der Waals forces. In 
addition, many extensions of the Standard Model that do not involve extra dimensions also 
predict the existence of new Yukawa or power-law forces. 

While gravity experiments at ranges > 10~ 3 m have found no convincing evidence of new 
forces or extra-dimensions, tests of Newtonian gravity over shorter separations were lacking 
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until recently. During the past few years a number of sub-millimeter gravity experiments 
were performed and stronger constraints on Yukawa corrections to Newtonian gravity for 



ranges ~ 10 4 m have been obtained 



HQ, HQ 



|. For significantly smaller separations, 



however, gravity loses its role as the dominant force acting between non-magnetic, electrically 
neutral interacting bodies. For these smaller separations, limits on new forces and extra- 
dimensions from force measurements must be extracted from the Casimir and van der Waals 
forces ^| which increase rapidly as the separation decreases. 

Improvements in the precision of Casimir force measurements coincided with the develop- 
ment of modern extra-dimensional theories. To date a number of Casimir force experiments 
have been performed using different techniques 19 



20 



21 



22 



23 



24, 
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26 , \2ji , and a 



precision of 2-3% (at the 95% confidence leve 



ration distances has 
e.g., in Refs. 0, 2U 



Deen achieved (see Ref. 



of the measured force at the shortest sepa- 
28]). We note that the 1% precision quoted, 



23[ corresponds to a 68% confidence level. To obtain good agreement 



between theory and experiment, it has been necessary to take into account corrections to 
the Casimir force due to finite conductivity of the boundary metals, surface roughness, and 



: to Unite conductivity oi t 

0, 0, 0, 



3. 



nonzero temperature 

Although these new precision measurements of the Casimir force were not especially de- 
signed or optimized to obtain stronger constraints on the predictions of extra-dimensional 
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38, 



This resulted in 



physics, some of them were used for this purpose |3j 
a strengthening of previously known constraints by a factor as large as 4500 in some re- 
gions within the range 10 _9 m to 10~ 4 m. This means that Casimir force measurements 
have become a powerful competitor to both accelerator and gravitational experiments in 
constraining theoretical models of high energy physics. In contrast to astrophysical and cos- 
mological constraints, the results of Casimir force measurements are less model dependent, 
reproducible, and therefore more reliable. 

In this paper we present detailed results of new, increased precision, Casimir force mea- 
surements between a Cu-coated plate and an Au-coated sphere (see Ref. 0] for preliminary 
data). The use of a microelectromechanical torsional oscillator (MTO) and of interfero- 
metric measurements of the sphere-plate separations permitted much higher sensitivity to 
be achieved than in previous Casimir force experiments. A careful error analysis has been 
performed, and the experimental precision was determined at the 95% confidence level. The 
complete theory of the Casimir force, taking into account finite conductivity and surface 
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roughness corrections has been applied to the experimental configuration. The finite con- 
ductivity corrections were computed by the use of the Lifshitz formula j^lj and tabulated 
optical data for the complex index of refraction. The surface roughness was modelled using 
atomic force microscope (AFM) images of the interacting surfaces. A comparison of the 
complete theory with experimental data shows that they are in agreement over the whole 
measurement range. A minor disagreement at the shortest separations noted in Ref. j3] is 
explained by the incomplete theory of roughness corrections used in the earlier analysis. 

Our results were sufficiently precise to shed light on the temperature-dependence of the 
Casimir force. Several theories predicting large thermal corrections to the Casimir force at 
sma „ sepals QQQ we,, elated for o„ r expert. It is we.. k aow„ that these 
corrections disagree with the results obtained for ideal metals in the framework of quantum 
field theory at nonzero temperature in the Matsubara formulation. Our experimental data 
support the results of Refs. j^, 0|, which are consistent with the conclusions of thermal 
quantum field theory, while ruling out the existence of large thermal corrections at small 
separations as proposed in Refs. [3,13 and in Ref. j^. 

Finally, the experimental results presented here are used to constrain the predictions 
of extra-dimensional physics in the nanometer separation range. The contributions from 
a Yukawa-type hypothetical force have been calculated for our experimental configuration, 
taking into account the effects of surface roughness. The agreement between theoretical and 
measured values of the Casimir force leads to the strengthening of the known constraints on 
Yukawa forces by a factor of up to 11 within the 56-330 nm interaction range. In contrast 
to some previous constraints derived from Casimir force measurements, it is possible here 
to quantify the confidence level of the obtained results. 

This paper is organized as follows. In Sec. II the experimental configurations used for both 
static and dynamic measurements of the Casimir force are described. Sec. Ill presents the 
experimental results with a discussion of their precision. In Sec. IV we calculate the Casimir 
force taking into account all relevant corrections. Sec. V is devoted to the determination 
of the theoretical precision, the comparison of experimental results with theory, and to the 
evaluation of alternative methods for taking into account thermal corrections. In Sec. VI 
we use our results to obtain constraints on hypothetical forces predicted by models of extra- 
dimensional physics and extensions of the Standard Model. We conclude with Sec. VII which 
summarizes all of our results. 



5 



II. EXPERIMENTAL ARRANGEMENT FOR STATIC AND DYNAMIC MEA- 
SUREMENTS 



In our experiment, the Casimir force between two dissimilar metals (gold and copper) 
was measured using a MTO operating in both static and dynamic modes. In the static 
regime the Casimir force between an Au-coated sphere and a Cu-coated plate of the MTO 
was measured. In the dynamic regime the vertical separation between the sphere and the 
plate was changed harmonically with time. This leads to a measurement of the ^-derivative 
of the Casimir force, which is equivalent to measuring the Casimir force per unit area, 
or the Casimir pressure, for a configuration of two parallel plates (see below). Note that 
the dynamic measurement technique is used here to measure the usual (static) Casimir 
effect. Hence, this measurement is unrelated to the so-called dynamic Casimir effect which 
arises from the velocity dependent Casimir forces or the creation of photons by the rapidly 
oscillating plates 3^|. 

When using mechanical oscillators to measure forces, one has to confront the coupling of 
the oscillator with environmental vibrations. Compared with cantilever oscillators, torsional 
oscillators are less sensitive to mechanical vibrations that induce a motion of the center of 
mass. Furthermore, the miniaturization of the oscillators yields an improvement in its quality 



factor and sensitivity 



25 



45j . It is consequently advantageous to use an electromechanical 



torsional oscillator to measure the Casimir force between two metals. 

The experimental arrangement is shown schematically in Fig. ^ The MTO is made of 
a 3.5 /im thick, 500 x 500 /im 2 heavily doped polysilicon plate suspended at two opposite 
points by serpentine springs, as shown in the inset of Fig. [21 The springs are anchored to a 
silicon nitride (SiN^) covered Si platform. When no torques are applied the plate is separated 
from the platform by a gap ~ 2 /im. Two independently contacted polysilicon electrodes 
located under the plate are used to measure the capacitance between the electrodes and the 
plate. For the MTO employed in this experiment, we calculated the torsion coefficient k = 
(wt 3 E S i/6L serp ) ~ 9.5 x 10~ 10 Nm/rad |4fij . where w = 2 /im is the width of the serpentine, 
t = 2 /im is its thickness, L serp = 500 /im its length, and E Si = 180 GPa is Young's modulus 
for Si. This value is in good agreement with the measured value k = 8.6 x 10~ 10 Nm/rad. 
The edges of the plate are coated with 1 nm of Cr followed by 200 nm of Cu. This layer of 
Cu constitutes one of the metals used in the measurement of the Casimir force. 
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The remainder of the assembly, as shown in Fig.^ consists of an Au-coated AI2O3 sphere 
that can be brought in close proximity to the Cu-coated plate. AI2O3 spheres with nominal 
diameters ranging from 100 /iin to 600 fim were coated with a ~ 1 nm layer of Cr followed by 
a « 203 nm layer of gold. The coated sphere used in the experiment was subsequently glued 
with conductive epoxy to the side of an Au-coated optical fiber, establishing an electrical 
connection between them. The sphericity of the AI2O3 balls, as measured on a scanning 
electron microscope (SEM), was found to be within the specifications of the manufacturer. 
As an example, a 600 /im-diameter ball was found to have an ellipsoidal shape with major 
and minor semi-axes of (298 ± 2) /iin and (294.0 ± 0.5) /im, respectively. For the sake of 
clarity we will refer to the ball as a sphere in the remainder of the paper. Deposition induced 
asymmetries were found to be smaller than 10 nm, the resolution of the SEM. The entire 
setup (MTO and fiber-sphere) was rigidly mounted into a can, where a pressure < 10 -4 ton- 
was maintained. The can has built-in magnetic damping vibration isolation and was, in 
turn, mounted onto an air table. This combination of vibration isolation systems yielded 
base vibrations with Az rms < 0.05 nm for frequencies above 100 Hz. 

The fiber-sphere assembly was moved vertically by the combination of a micrometer- 
driven and a piezo-driven stage. The MTO was mounted on a piezoelectric driven xyz stage 
which, in turn, is mounted on a micrometer controlled xy stage. This combination allows 
positioning the Au-coated sphere over the Cu-coated plate. The separation between the 
sphere and the Si-platform was controlled by the z-axis of the xyz stage. A two color fiber 
interferometer-based closed-loop system was used to keep Zi constant. The error in the 
interferometric measurements was found to be Az^ ms = 0.32 nm, dominated by the overall 
stability of the closed-loop feedback system. Since this error is much greater than the actual 
mechanical vibrations of the system, the closed loop was turned off while data acquisition 
was in progress. 

The separation z meta i between the two metallic surfaces (see Fig. [l]) is given by z me tai = 
Zi — zq — Zg — bd, where b is the lever arm between the sphere and the axis of the MTO, and 
9 is the angle between the platform and the plate {6 <C 1 has been used), zq is the distance 
the bottom of the sphere protrudes from the end of the cleaved fiber, and z g includes the 
gap between the platform and the plate, the thickness of the plate, and the thickness of the 
Cu layer. An initial characterization of z , by alternately gently touching the platform with 
the sphere and the bare fiber yielded z = (55.07 ±0.07) /im. Also, z g = (5.73 ± 0.08) /im, 
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was determined using an AFM. Since errors in z g and z propagate to z meta i, it is necessary 
to provide a better characterization, which is described below. 

A force F(z) acting between the sphere and the plate produces a torque r = bF(z) = k9 
on the plate. (In all cases reported in this paper 9 < 10~ 5 rad, so C 1.) Under these 
circumstances 9 oc AC = C rig ht — Ci e f t , where C rig ht ifluft ) is the capacitance between 
the right (left) electrode and the plate (Fig. [T]). Consequently the force between the two 
metallic surfaces separated by a distance z is F(z) = kAC, where k is the proportionality 
constant. The capacitance was measured using the circuit schematically shown in Fig. El 
The 10 nV/Hz 1//2 electronic noise of the amplifier stage is equivalent to an angular 
deviation 59 ~ 1(T 9 rad/Hz 1 / 2 . This noise level is much smaller than the thermodynamic 



3, 



s, 



1/2 



3 x l(T 7 rad/Hz 1/2 , (1) 



noise at the measuring frequency / = 10 kHz 3> fo = 687.23 Hz 

lk B TQ //o\ 2 l V2 
t/o« V / / 

where ks is Boltzmann's constant, T = 300 K, and the quality factor Q ~ 8000. (See the 
inset of Fig. [TJ) Using these values a force sensitivity 5F = kS X J 2 jh ~ l^pN/Hz 1 / 2 is 
obtained. 

To provide a calibration of the proportionality constant k between F and AC, we applied 
a known potential difference between the Au-coated sphere and the Cu film. This was done 
at z m etai > 3 jum, in configurations where the contribution from the Casimir force is smaller 
than 0.1 % of the total force. Thus, the net force can be approximated by retaining only 



the electrostatic force F e between a sphere and an infinite plane 

F e = 2ne (V Au -V o r± ^^1 . (2) 

L — ' smh(nu) 

n=l v ' 

Here eo is the permittivity of free space, Va u is the voltage applied to the sphere, Vo is 
the residual potential difference between the metallic layers when they are both grounded, 
and coshu = [1 + z/R], with z = z meta i + 2<5 (R is the radius of the sphere). 25 is 
the average separation between the metal layers when the test bodies come in contact (so 
that z meta i = 0), and is primarily determined by the roughness of the films. Although the 
profile of the roughness does not affect the value of the electric force at large separations, 
its presence should be taken into account in the determination of the separation between 
the smoothed out surfaces. In Eq. (J2J) it was found that only the first two terms of the z/R 
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expansion gave a significant contribution. Fig. |3] shows the dependence of 9 (and hence the 
electrostatic force) on the applied voltage Va u - The minimum in the force was found to 
occur when Vo = (632.5 ± 0.3) mV, which reflects the difference in the work functions of the 
Au and Cu layers. This value was observed to be constant for z in the 0.2-5 /zm range and 
it did not vary when measured over different locations in the Cu layer. 

Once the value Vo was found, Eq. (j2J) was used to determine three different parameters: 
(i) the proportionality constant k between the sphere-plate force and the measured difference 
in capacitances; (ii) the radius R of the coated sphere; and (iii) the increase in the separation 
2<5o between the two metallic layers. Two of the curves obtained in the 3-5 /zm range are 
shown in Fig.0] A simultaneous fit to more than 100 such curves yields k = (50280 ±6) N/F, 
R = (294.3 ± 0.1) /mi, and 5 = (39.4 ± 0.3) nm. 

The above force sensitivity can be improved by performing a dynamic measurement which 



directly uses the high quality factor of the MTO |25j]. In this approach, the separation 
between the sphere and the oscillator was varied as Az meta i = Acos(cu r t), where io r is the 
resonant angular frequency of the MTO, and A was adjusted between 3 and 35 nm for values 
of z me tai between 0.2 and 1.2 /zm, respectively. The solution for the oscillatory motion yields 

. . r h 2 fiF„~\ 

(3) 



2 1 



u ore 
Iuq dz 



where ujq ~ y/~K/I for Q ^> 1, I ~ 4.6 x 10 _17 kg-m 2 is the moment of inertia of the 
oscillator, and Fq is the Casimir force between the sphere and the plate. Since A <C z meta i, 
terms of higher order in dFc/dz introduce a ~ 0.1% error at separations z > 250 nm. As 
before, Eq. (J2J) was used to calibrate all constants. We found ujq = 2ir x 687.23 Hz, and 
b 2 / I = 1.2978 x 10 9 kg _1 . With an integration time of 10 s using a phase-lock-loop circuit 
50j, changes in the resonant frequency of 10 mHz were detectable. 

The main source of error for the frequency measurement, 5f = 10 mHz, is the error in the 
trigger of the frequency meter. This error 5t originated in the jitter of the signal due to the 
relatively large thermodynamic noise- induced 59 observed at resonance 48] . This dominant 
noise is 



V=^* (4) 



where T is the integration time. 
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Unlike the static regime where forces are measured, in the dynamic regime the force 
gradient dFc / dz is measured using Eq. (J3J) by observing the change in the resonant frequency 
as the sphere-plate separation changes. According to the proximity force theorem j5ll l52|. 

F c (z) = 2nRE c (z), (5) 

where Ec(z) is the Casimir energy per unit area for two infinitely large parallel plates com- 
posed of the same materials as the sphere and plate (see Sec. IV for details). Differentiating 
Eq. (jSJ) with respect to z one obtains 

where Pc{ z ) is the force per unit area between two infinite plates. Thus, in the dynamic 
regime the Casimir force gradient between the sphere and plate can be directly related to the 
Casimir pressure between two infinite parallel plates composed of the same materials. Since 
the formula for the Casimir pressure for the parallel plate configuration is readily obtained, it 
proves more convenient to use Eq. (jHJ) to express the results for the dynamic measurements 
in terms of the parallel plate pressure Pc instead of the force derivative dFc/dz. From 
Eq. Q , the error in the measurement of the frequency translates into an equivalent pressure 
sensitivity given by 5P C ~ 4 x 1CT 4 Pa/Hz 1//2 . 

Finally, one additional test was performed, in this case to analyze the influence of the 
finite extent of the Cu layer on the measured forces. The analysis was done using, once again, 
the electrostatic force given by Eq. (J2J). Fig.Elshows the relevant data. When the sphere was 
moved parallel to the axis of the oscillator over 20 /mi no change in 9 was observed. When 
the motion was instead perpendicular to the axis of the oscillator the dependence expected 
from Eq. (J2J) was obtained within the experimental error. We thus conclude that there are 
no significant deviations from the assumption that the Cu plane is of infinite extent. 



III. EXPERIMENTAL RESULTS AND THEIR PRECISION 



A. Surface roughness of samples 

The electric force measurements used to calibrate the apparatus were performed at large 
separations so that the Casimir force could be neglected, and the plate-sphere separation was 
calibrated between the middle levels of the surface roughness. By contrast, the Casimir force 
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at small separations depends sensitively on the profile of the surface roughness, and hence 
the surface roughness should be carefully analyzed and characterized. The topography of 
the metallic films was investigated using an AFM probe with a radius of curvature r c = 5 nm 
in tapping mode. Regions of the metal plate and the sphere varying in size from 1 /zm x 1 /j,m 
to 10 /iin x 10 fim were scanned. A typical surface scan of a 1 fim x 1 /xm region is shown in 
Fig. where the lighter tone corresponds to higher regions. As seen in Fig. EJ the major 
distortions are the large mounds situated irregularly on the surface. Also noticeable in Fig. El 
are streaks which arise from high frequency noise with amplitude h rms ~ 1 nm. This noise 
is caused by the oscillation of the free-standing MTO while acquiring AFM images. 

In order to include the effects of surface roughness in the Casimir force calculations, 
the fraction of the surface area Vi with height hi is needed. Data resulting from the most 
representative scan of a 10 /j,m x 10 /im region are shown in Fig. [7| The heights hi are plotted 
along the vertical function of the fraction Wi of the total surface area having height 

less than hi + \. The width of each horizontal step is equal to the fraction of the total 
area Vi with heights hi < h < hi+\. For example, regions with heights h < h2 = 8.2 nm 
are referred to the first (bottom) distortion level hi = with v\ = 0.0095, and regions 
with heights h 2 < h < h 3 = 9.87 nm are referred to the second distortion level h 2 with 
t>2 = 0.0121. Regions with heights h > h§g = 98.5 nm are referred to /159, and occupy a 
fraction t> 59 = 0.085 of the total area. Evidently tOj = v 1 + v 2 + ■ ■ ■ + v j, and w 59 = 1. 

The data of Fig. will be used in Sec. IV for the computation of the Casimir force taking 
into account roughness corrections (in Ref. j3] a simplified model of roughness was used). 
These data are also required for the precise determination of the so-called zero roughness 
level Hq relative to which the mean value of the function, describing the total roughness, 
is zero. For convenience in comparison with theory, all separation distances in the Casimir 
force measurements presented below are measured between the zero roughness levels. 

The zero roughness level Hq is found using the equality 



By combining the data of Fig. 0with Eq. (|7j) . one finds Hq ~ 35.46 nm. It is seen that the 
zero roughness level is slightly different from the quantity 5q ~ 39.4 nm, the correction to 
the separation on contact determined by the results of the electric force measurements (see 
Sec. II). The separation distances between zero roughness levels z — d — 2H used below (d 



59 




(7) 



i=l 
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is the separation between the bottom roughness levels) are larger by 2(5q — Hq) = 7.88 nm 
than the separations d — 25q defined by the electric force measurements. (The difference 
between H and Sq can be explained by a minor modification of the highest roughness peaks 
before the beginning of the Casimir force measurements.) 

B. Static measurements of the Casimir force 

To measure the Casimir force in the static regime, the bridge (schematically shown in 
Fig. |2I) was first balanced using two identical capacitors replacing the MTO. Then the MTO 
was put back in place and the voltages V\ and V2 adjusted to give a null signal. This last 
adjustment is required to take into account residual asymmetries in the MTO. It was found 
that the difference between V\ and V2 corresponded to a variation in z, 5z < 6nm. Once 
the MTO was mounted and the can was evacuated, the sphere was brought into proximity 
with the MTO and the electrostatic measurements were performed. Without breaking the 
vacuum in the system the Casimir force measurements were then carried out. Fig. |H1 shows 
one such data set out of 19 runs. Each data point was obtained with an integration time of 
10 s, a time interval which represents a good compromise for the ~ 300 data points taken 
per run. It is worth mentioning, however, that the force sensitivity can be improved by 
using longer integration times. 

The following analysis of the experimental precision is based on all n — 19 series of 
measurements. For this purpose we calculate the mean values of the measured Casimir force 

1 n 

rw = -E^w ( 8 ) 

i=l 

at different separations z m within the measurement separation range from 190 nm to 
~ 1.15 /zm. The mean square error of F^ xp is equal to 

"£[F^(z m )-F°J(z m )} 2 \ . (9) 
i=i J 

Our calculations show that s n (z m ) does not depend sensitively on m. The largest value 
s n = 0.143pN can be taken as the value of the mean square error of F^ xp within the whole 
measurement range. Taking into account that at a = 95% confidence level the Student's 
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coefficient is t a>n = 2.1, one obtains for the random absolute error of the Casimir force 
measurements in the static regime 

A W F ex P = ^ w ^ 03pN _ (1Q) 

In fact this effectively gives the total absolute error ££°IFq XP since the systematic errors are 
far below 0.1 pN for an integration time of 10 s. As a result, at the shortest separation of 
about 188 nm the relative error of the Casimir force measurement is 0.27%. We note that 
the true value of the Casimir force at a separation z lies in the confidence interval 

with a probability of 95%. 

C. Dynamic measurements of the Casimir pressure 

The results for the parallel plate Casimir pressure in the dynamic regime are shown in 
Fig. El which presents one data set of 5 runs. Although the measurement was extended 
down to a separation of ~ 180 nm, only data for separations above 260 nm are plotted. Data 
points below 260 nm show effects of nonlinear behavior of the oscillator 2^| . 

The error analysis is performed in the same way as for the static regime. For I = 5 runs, 
with m 300 points per run, the largest mean square error is found to be si = 0.11 mPa. 
The Student's coefficient at a = 95% confidence level is t a j = 2.8 leading to a random 
absolute error ^^P^ xp = 0.31 mPa. In the case of dynamic measurements, however, it is 
not possible to neglect systematic errors compared to random errors. As noted in Sec. II, 
there is an error 5u ~ 2tt x 10~ 2 Hz which from Eq. leads to the absolute error in 
dFc/dz « 4.2 x 10~ 7 N/m. The latter, when combined with the error associated with the 
sphere's radius (equal to 0.1 ^m), leads via Eq. (JHJ) to a systematic error on the pressure 
^dyn^ > (j V i w hich varies from 0.31 mPa at the shortest separation z = 260 nm to 0.23 mPa at 
all separations z > 450 nm. As a result, the total absolute error of pressure measurements 
in the dynamic regime is equal to 

xtot r> ex P _ Arand r>exp , \syst r>exp / 1 1 \ 

^■dyn r C — ^dyn r C + ^dyn^C ■ \ lL ) 
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This error is ^-dependent and varies from 0.62 mPa at z = 260 nm to 0.54 mPa at z > 450 nm. 
Hence, the relative error of the Casimir pressure measurement at the shortest separation of 
about 260 nm is 0.26%. 



IV. THEORETICAL DETERMINATION OF THE CASIMIR FORCE 
A. Casimir force and pressure including finite conductivity 

As described in Section II, the static and dynamic measurements were carried out using a 
Au-coated sphere over a Cu-coated plate. The thicknesses of both metal coatings were much 
larger than the plasma wavelength X p of both metals so that we can calculate the Casimir 
force as if the sphere and plate were composed of solid Au and Cu respectively. As a first 
approximation, we consider the plate to have an infinite area. (Corrections due to the finite 
size of the plate will be estimated below.) The Casimir force between an infinite plate and 
a sphere (which was measured in the static regime, see Sec. II) can then be found using the 

nn 

Lifshitz formula for two plane parallel plates and the proximity force theorem |41l . I53IJ 

p poo POO 

fcW = -| k A J d^{\n[l-r^^k ± )rf\^k ± )e~ 2 ^\ (12) 

+ In [l - r^fo fc ± )ri 2) (£, k^e^] } . 

Here q 2 = k\ + £ 2 , k± is the modulus of the wave vector in the plane of the plates, and 
rji Z j_ (I = 1, 2 for Cu, Au, respectively) are the reflection coefficients for two independent 
polarization states computed along the imaginary frequency axis u = iC,. We note that the 
errors introduced by the proximity force theorem used to derive Eq. (fT2l) are smaller than 

an 
|5^. This is a correction of less than 0.06% at the shortest separation distance, 

where the force measurements are most precise. (Recall that the experimental precision at 

z = 188 nm was found to be 0.27%.) Eq. ()12|) takes into account the finite conductivity 

corrections to the Casimir force due to real metal boundaries, but does not consider the 

effect of surface roughness and thermal corrections which will be treated later. 

The reflection coefficients used in Eq. ()12|) can be represented in terms of either the 

dielectric permittivity or the surface impedance. In terms of the dielectric permittivity, as 

in the original Lifshitz formula (denoted by the subscript L), the reflection coefficients are 
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given by 



feCO 



r(0| 



*0e 



.(0 



9 



_jfe(0 



fc(0+e(0(i£)g' '^^''^ g + jfe(0 ' 
where fc^ 2 = k\ + £^(i£)£ 2 - I 11 terms of the surface impedance, the reflection coefficients 
are given by |34j 



(13) 



z (t) (*flg-g 



(14) 



where Z^ l \i^) is the impedance 
It has been shown recently 



'or Cu or Au computed along the imaginary frequency axis. 
| that for real metals at nonzero temperature and interacting 
at relatively large separations, the representation described in Eq. (ftl^) is preferable. It is free 
of contradictions with thermodynamics, which arise when reflection coefficients expressed 



by Eq. (|13|) are used in combination with the Drude dielectric function containing a nonzero 



relaxation parameter 56|. Within the experimental separation ranges of Sec. II (0.2-1.2 
/zm) the characteristic angular frequency of the Casimir effect u c — 1/ (2z) lies within the 
region of infrared optics (using c = 1), where the dielectric permittivity of the plasma model 
and the surface impedance along the imaginary frequency axis are given by 



to r 



|57| , both formulations of the reflection coefficients, Eqs. (fT3|) and ([I 



(15) 



As demonstrated in 

lead to the same computational results for the Casimir force, Eq. (f*H?|) . when Eq. (|15|) is 
taken into account. 

In the dynamic regime of Sec. II, the Casimir force gradient between the sphere and plate 
was measured. As shown above, this force gradient can be re-expressed in terms of the 
Casimir pressure acting between two parallel plates composed of the same materials. This 
Casimir pressure between parallel plates can also be represented in terms of the reflection 
coefficients 



Pc(z) 




- 1 



(16) 
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For separations z > X p the Casimir pressure can be computed using either set of reflection 
coefficients, Eq. (fT^j) or Eq. (JHj), with or Z(i£) given by Eq. (JT3J). At the shortest 

separations z < X p the impedance at characteristic frequencies is not small, and Eq. (J13j) 
should be used to calculate both the Casimir pressure and force. The most accurate results 
at these separations are obtained by using tabulated data for the imaginary part of the 
dielectric permittivity [58] . These data are substituted into the dispersion relation 



2 f°° ojlmeju ) 

to obtain the dielectric permittivity along the imaginary frequency axis. The Casimir force, 

L This procedure 



Eq. (jl2)l . and pressure, Eq. ()16j) . can then be calculated as in Refs. Q,l2 
was used here to calculate Fc{z) and Pc(z). The available tabulated data 



were extended 



using the Drude model with the following plasma frequencies and relaxation parameters 



(1) - 9.05 eV, 4 2) = 9.0 eV, 7W = 30meV, 7^ = 35meV, for Cu and Au, respectively. At 



large separations z > X p our results almost coincide with those obtained in the framework 
of the plasma model given by Eq. (fTKJ) . Note that at smaller separations z < X p our results 
are practically independent of the chosen extrapolation and are completely determined by 
the available tabulated data. 

It is useful to compare the calculated results for the real metals used with the ideal case 
when both sphere and plate were composed of a perfect metal with infinite conductivity. For 
example, the ratios of the calculated Casimir force to the force between an ideal sphere and 
plate are 0.467, 0.544, and 0.842 at separations z = 70 nm, 100 nm, and 500 nm, respectively. 
Similarly, the ratios of the calculated Casimir pressure between real plates to the pressure 



between ideal plates are 0.393, 0.468, and 0.799 for t 
very close to those computed for Au-Au and Cu-Cu 
similarities of the optical data for Cu and Au. 



re same separations. These results are 
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which can be explained by the 



B. Casimir pressure and force including surface roughness 

As mentioned previously, the theoretical results obtained from Eqs. ()12|) and ()16|) take 
into account the finite conductivity of the boundary metals but neglect the effect of surface 
roughness. This effect, however, may constitute a correction of several tens of percent at 
the shortest separations depending on the character of the roughness. Hence, a precise 
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computation of the Casimir force requires a careful characterization of the roughness as was 
performed in Sec. IIIA. 

The roughness correction to the Casimir interaction is computed using the AFM images of 
the surfaces, like the one shown in Fig. El As seen from Fig. |Hl the characteristic longitudinal 
scale of the surface roughness is larger than the surface separation (especially in the region 
of the shortest separations where the effect of roughness is most significant). In this case 
the additive method j^, |3| can be used to calculate the Casimir force taking account of 
roughness. As a result, the Casimir pressure between two plates with roughness corrections 
taken into account is given by 



where Vi is defined in Sec. IIIA, Pc( z ) is given by Eq. (flBjl. and the index i relates to one 
plate and j to the other. A plot of the relief heights hi (i = 1, 2, . . . , n = 59) on the plate 
versus the fraction of the plate area Wi with height h < hi + \ is shown in Fig. It should 
be remembered that H = 35.46 nm is the zero roughness level from which all separations z 
are measured (see Sec. IIIA). Note also that the separation distances z + 2H — hi — hj in 
Eq. (|TH|) may be much smaller than X p for large relief heights hi, hj. Hence, one should use 
the optical tabulated data when calculating the finite conductivity corrections. 

In a similar manner, the Casimir energy for two parallel plates with roughness is given 
by Eq. (fTHj) after substituting Ec for Pc- It then follows that by the use of the proximity 
force theorem, the Casimir force for a sphere above a plate, accounting for both roughness 
and finite conductivity corrections, is given by 



where F c (z) is given by Eq. (fl"2"|) . 

We note that Eqs. (j!8|) and (|T9"|) describe not only the separate effects of finite conductivity 
and surface roughness on the Casimir pressure and force, but also their combined effect. This 
is especially important at the shortest separations where the corrections are not small and 
cannot be represented as a product of two separate factors, one each for roughness and finite 
conductivity. 



n 




(18) 



n 




(19) 
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V. COMPARISON OF THEORY WITH EXPERIMENT AND TESTS OF ALTER- 
NATIVE THERMAL CORRECTIONS 



The theoretical Casimir force Fc(z) acting between a sphere and a plate was computed 
using Eq. (|19|) for all separations where it was measured (19 sets of measurements containing 
~ 300 experimental points each). In Fig. [IJ3 the difference between the theoretical and 
experimental force values 



as a function of surface separation Zj is presented for one set of measurements. As can be 
seen from the figure, the values of AFc{zi) are clustered around AFc(zi) = 0, demonstrating 
good agreement between theory and experiment. 

To quantify the level of agreement between theory and experiment we consider the root 
mean square (r.m.s.) deviation afj defined as 



where N is a number of points under consideration. 

For example, if the first 250 points from all 19 sets of static measurements are considered 
(separations larger than 1 /jm are not considered due to the large experimental relative error), 
one obtains N = 4750, and cr^ 50 pa 0.6 pN, which is less than the theoretical error (see below) 
but two times larger than the absolute error of the force measurements given by Eq. (jlOj) . 
The r.m.s. deviation depends slightly on the separation region under consideration. Thus, 
if separations z > 400 nm are considered, the first 185 points of all sets of measurements 
lead to N = 3515, cr][ 515 pa 0.45 pN. For one set of measurements shown in Fig. ITUl iV = 250 
and <7<^ pa 0.66 pN (close to the above value when all sets of measurements are considered). 
It can also be seen from Fig. E3 that there is a slight shift of the mean difference force value 
below the zero level equal to — 0.045 pN, whose magnitude is much smaller than the absolute 
error of the force measurement. 

A similar analysis was performed for the parallel plate Casimir pressure obtained from the 
dynamic measurement. The difference between the theoretical and experimental pressures 
APc(zi) is defined as in Eq. (J2UJ) . and the r.m.s. deviation is defined as in Eq. (l2Tj) 




(20) 
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with the substitution of Fc for Pc- Considering the first 235 points of all 5 sets of dynamic 
measurements one obtains N = 1175, and cr^ 75 w 0.5 mPa. Once again, the r.m.s. deviation 
depends slightly on the separation interval. Thus, for 310 nm < z < 420 nm (35 points), 
N = 175 and a± 75 ~ 0.44 mPa if all five sets of measurements are considered. For z > 310 nm 
(220 x 5 = 1100 points), of 100 « 0.34 mPa. 

In Fig. the quantity APo(zi) is presented for one set of measurements (N = 235). 
Here the shift of the mean difference pressure value below zero is equal to —0.26 mPa (whose 
magnitude is less than the absolute error of pressure measurements). The r.m.s. deviation 
between theory and experiment in Fig. ^^is cr^ 5 rs 0.43 mPa, i.e., even less than for all five 
sets of dynamic measurements. Both Figs. El and ^2 demonstrate good agreement between 
theory and experiment. 

To estimate the theoretical precision, we consider other effects, in addition to finite con- 
ductivity and surface roughness, which were not taken into account in Eqs. (fTHjl. (fT9")l for the 
Casimir pressure and force. Most prominent among these are finite-temperature corrections, 
which will be dealt with in detail later in this section. Another correction which should be 
considered arises from the fact that the plate used in the experiment is not infinite in ex- 
tent. For the plate used, whose dimensions are 500 x 500 /mi 2 , this correction is evidently 
negligible if the sphere's center is located above the mid-point of the plate. However, in 
the configuration shown in Fig. the sphere's center is above a point displaced from the 
right plate boundary by a distance L = 50 /an. In this case, the Casimir force should be 
multiplied by a correction factor (3 that was found in Ref. [6jj 



where R is the radius of the sphere. In our case, at the largest separation z = 1 /im, this 
leads to (3 = 0.9865. However, taking into account the fact that only one of the four sides of 
the plate is close to the sphere center projection reduces the correction factor to (3 = 0.997. 
Thus, the correction is in fact less than 0.3%, and at a separation z = 500 nm it is less 
than 0.04%. The same considerations apply to the pressure between two parallel plates. 
Consequently, the correction to the theoretical force and pressure due to the finiteness of 
the plate is much smaller than the uncertainty introduced by the sample-to-sample variation 
of the optical data. These sample-dependent variations in the index of refraction may lead 




(22) 
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to an error of approximately 1% 29, 30]. An even smaller uncertainty is introduced into 



the computation by the effect of roughness when the topography of the surface is carefully 

n. 

characterized. A recently discussed correction due to the surface plasmon [61] is also not 
significant. The surface plasmon propagates when the frequency of the electromagnetic wave 
is greater than uo p . Recall that in our case the shortest separations are 190 nm and 260 nm 
in the static and dynamic regimes, respectively. This leads to the highest characteristic 
frequencies uj c = 0.8 x 10 15 rad/s and 0.6 x 10 15 rad/s, respectively, which are 17 and 23 
times, respectively, less than the plasma frequency for Au and Cu. As a result, the correction 
due to the surface plasmon in our case is much less than 1% even at shortest separations 
used. 

As mentioned earlier, an important correction which may influence the magnitude of the 
Casimir interaction is due to the effect of nonzero temperature. This has been the subject 
of —able c_ y during the U-t few years (see, e,„ Q Q Q Q Q QQ 



0,0). For an ideal metal the thermal correction can 
formulation of thermal quantum field theory (TQFT) [31 



)e determined using the Matsubara 
0]. The question is 



0, 0,00, 



whether the thermal correction for good (but real) conductors is small at small separations, 



as is the case for idea 
The approach of Refs. 



metals, or does it differ qualitatively from that for ideal metals. 
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which we call "traditional" since it yields results consistent 
with earlier studies of thermal effects, leads to a qualitatively identical thermal correction 
for real and ideal metals, as given by TQFT. For example, in the configuration of two 
parallel plates made of Cu and Au, the Casimir pressures at z — 300 nm and z = 500 nm 
at zero temperature are equal to — 136 mPa and — 17.0mPa, respectively. The traditional 



thermal corrections 
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331 ] at room temperature T = 300 K for these cases are equal to 



—0.00863 mPa and —0.00441 mPa, respectively. (For comparison, in the case of ideal metals, 
the corresponding Casimir pressures at T = OK are equal to — 160.43mPa and — 20.79 mPa, 
respectively, and the thermal correction at T = 300 K is — 0.00204 mPa and is independent of 
separation.) It is clearly seen that the traditional thermal corrections at room temperature 
are very small, and even the improved sensitivity of our experiment is not sufficient to 
measure them. Their contribution to the Casimir pressure is ~ 0.006% at z — 300 nm and 
~ 0.03% at z = 500 nm, and can therefore be neglected. 



The situation with the alternative thermal corrections, proposed in Refs. 



42 



Q,Q, 



is quite different. These corrections at separations z > X p are much greater than those 
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predicted by the traditional approach. According to Refs. 3, Q , in the case of two parallel 
plates made of real metals there is an additional thermal correction linear in temperature 
given by 



where 



A T Pg\z) 



k B T 
16nz 3 



V 2 dy 



(2), 



- 1 



(23) 



ri°(0,y) 



y — V y 2 + ^z 2 u : 



y+ \/y 2 + ^z 2 uj { i y 



(24) 



The effect of this correction is not as small as for the traditional approach. For example, 



at separations z = 300 nm and 500 nm Eq. (|231) leads to A r p£ X) = 4.89 mPa and A T P C - 
1.23 mPa, respectively, at T = 300 K, i.e., to 3.6% and 7.24%, respectively, of the total 
Casimir pressure. Since these values far exceed the errors of the present experiment, the 



(i) 



new resu 
in Refs. 



ts can be used as a decisive experimental test for the theoretical predictions made 
42IU. 



In Fig. El the difference between the theoretical and experimental Casimir pressures 



pth,i _ pexp - g p resen t ec i as a function of separation for the same set of measurements as used 



in Fig. ^2 Here, however, the theoretical pressure Pq' 1 is computed using the alternative 
thermal correction given by Eq. ()23|) . It is obvious that at separations < 700 nm the quantity 



pth,i _ pexp (jgyjg^gg significantly from zero. At the shortest separation z = 260 nm this 
deviation reaches 5.5 mPa. Thus, the linear thermal correction to the Casimir pressure 
proposed in RefsJiJ Q is ruled out by the present experimental results. (Note that in 
recent preprint 68] qualitative arguments are presented on the role of the finite size of the 
plate and finite thickness of the gold layer. According to j3| these effects could lead to a 
25% discre pan cy between the experiment of Ref. |19'], and the proposed alternative thermal 
correction at a separation 1 /im.) 

A second alternative thermal correction to the Casimir pressure between real metals was 
proposed in Ref. j^. For our case of two metal plates it is expressed as 



A T Pg\z) 



C(3) + A T P«(z) 



(25) 
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where ((x) is the Riemann zeta-function, and AtP^\z) is defined in Eq. ([23)1 . This cor- 
rection at T = 300 K is also much larger than the traditional one for the separations used 
in our experiment. At the separation z = 300 nm Eq. ()25|) leads to AjP^ = — 2.44mPa, 
i.e., to a 1.8% correction to the total Casimir pressure (the experimental precision at this 
separation is 0.43%). Thus, the present experiment also provides a test for the theoretical 
predictions of Ref . |44| . 

In Fig. ED the difference between the theoretical and experimental Casimir pressures 
pth,2 _ pexp j g p resen ^ ec [ as a function of separation for the same set of measurements as 
in Fig. ^2 Unlike Fig. El the theoretical pressure P<^' 2 is computed using the second 
alternative thermal correction given by Eq. (|25|) . As seen from Fig. El at separations less 
than 600 nm the quantity P<^' 2 — P^ p deviates significantly from zero, and reaches 5mPa 
at a separation z = 260 nm. It follows that the thermal correction of Ref. Q| is also in 
contradiction with the results of the present experiment. 

The conclusion that the alternative approaches to thermal corrections to the Casimir 
force do not agree with our experimental results is not surprising upon recognizing that the 
approaches of Refs. (2, Q,Q| violate the Nernst heat theorem pfl |. To correctly describe 
the influence of thermal effects on the Casimir force between real (i.e., non-ideal) metals 
requires a proper understanding of the zero-frequency contribution in the Lifshitz formulas 
given by Eqs. (|T2*j) and ()16|). When the Drude dielectric function with a nonzero relaxation 
parameter is substituted into Eq. (jl3j) . the approach of Refs. , 43 ] follows and the thermal 
correction given by Eq. (}23|) is found. If one then modifies the reflection coefficient r i (0, k \ ) 
by setting it equal to unity (as for real photons), the thermal correction (j25j) is obtained [44]. 
However, to avoid contradictions with fundamental physical principles and experiment, it is 
necessary to start from the physical behavior of the surface impedance in the appropriate 
range of characteristic frequencies (i.e., infrared optics in our case), and to extrapolate to 
zero frequency . Under these conditions, the traditional thermal correction between real 
metals is recovered, and this thermal correction is in agreement with both thermodynamics 
and the present experiments. Furthermore, it transforms smoothly into the limiting case of 
ideal metals as described by the Matsubara formulation of TQFT. 

To conclude this section, the theoretical uncertainty, which is approximately 1% of the 
Casimir force or pressure, exceeds the experimental uncertainty at the shortest separations 
used in this experiment (0.27% for the force at 188 nm and 0.26% for the pressure at 260 nm). 
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However, at separations z > 370 nm the experimental uncertainties exceed the theoretical 
uncertainties. Within the achieved levels of precision there is good agreement between theory 
and experiment over the entire measurement range. 



VI. CONSTRAINTS ON NEW YUKAWA FORCES AND EXTRA- 
DIMENSIONAL PHYSICS 

As mentioned in the Introduction, in many extensions to the Standard model, including 
theories with large compact extra dimensions the potential energy between two point 
masses mi and m 2 separated by a distance r is given by the usual Newtonian potential with 
a Yukawa correction 



parated d; 



V(r) = -^p(l + ae-^), (26) 

where a is a dimensionless constant characterizing the strength of the Yukawa force, and A 
is its range. For theories with n > 1 extra dimensions a ~ 1 and A ~ R n , where R n is the 
size of the compact dimensions, Eq. (J2l)|) holds under the condition [j| 

2/n 



" >>R "^[W ~ 1 " m ' (2?) 

For n — 1 it follows from Eq. ()27|) that Ri ~ 10 15 cm which is excluded by solar system tests 
of Newtonian gravity If, however, n = 2 or n = 3, the sizes of extra dimensions are 
R 2 ~ 1 mm or R 3 ~ 5 nm, respectively. While recent gravity experiments have investigated 
millimeter distance scales without finding evidence of new physics, gravity remains poorly 
tested at scales < 10~ 4 m (see Ref. (ill for a review). 

For models of non-compact (but warped) extra dimensions [2| the potential energy takes 
the form of the Newtonian potential with a power-law correction 



™-J5S=Sfl + a k), (28, 



where r>l/^ and 1/k is the warping scale. The correction in Eq. (|2*Hj) can be generalized 
to arbitrary inverse powers, 



Gm 1 m 2 
Vi(r) = 



23 



1 + an — 
V r 



(29) 



where on is a dimensionless constant, I is a positive integer, and tq = 10~ 15 m. 

We note that Yukawa and power-law corrections given in Eqs. ()26|) and (}2*9"j) also arise in 
ways unrelated to extra-dimensional physics. For example, the Yukawa potential describes 
new forces generated by the exchange of light bosons of mass \i — 1/A, such as scalar axions, 
graviphotons, hyperphotons, dilatons, and moduli among others (see, e.g., 00,0,0, 



731]). For such forces the interaction constant a could be much larger than unity. Power-law 



corrections, as in Eq. (}2*9l. arise from the simultaneous exchange of two photons or two 
massless scalars (/ = 2 [rJ]). two massless pseudoscalars (/ = 3 jz^l, Q)' and from the 
exchange of a massless axion or a massless neutrino-antineutrino pair (1 = 5 jzf| 3]). 

The agreement between theory and experiment for our Casimir force measurements can be 
used to set new constraints on the Yukawa strength a as a function of A from Eq. (J26j) . The 
total force acting between a sphere and a plate due to the potential described by Eq. (|2T)j) 
is obtained by integration over the volumes of the sphere and the plate, and subsequent 
differentiation with respect to z. In fact, the contribution of the Newtonian gravitational 
force is very small and can be neglected. To prove this, let us consider a sphere above the 
center of an enlarged plate modelled by a disk with a radius Lq 3> R = 294.3 pm. (In the 
present experiment the projection of the sphere center is displaced by 450 /im from one edge 
of the plate, and by 50 /im from the other). In this case the Newtonian gravitational force 
is given by 39] 

8 f D R\ 

F N ~ --iT 2 Gp disk p sphere DR 3 1 - — - — ) , (30) 
6 \ ZLq Lq J 

where D = 3.5 /im is the thickness of the plate. To obtain an upper limit, let us neglect 
the layered structure of both test bodies and set L —> oo, pdisk = Pcu = 8.93 x 10 3 kg/m 3 , 
and Psphere = PAu = 19.28 x 10 3 kg/m 3 . It then follows from Eq. (jHSJ) that F N »s —3.2 x 
10~ 17 N. This value is four orders of magnitude smaller than the absolute error of the force 
measurement in the static regime. Hence, the contribution of the Newtonian gravitational 
force can be neglected at this stage. 

For a Yukawa force between a sphere and a plate, the constraints should be calculated 
considering the detailed structure of the sphere and plate. The sphere of density p s = 4.1 x 
10 3 kg/m 3 was coated with a layer of Cr of thickness Acv = 1 nm with po = 7.19x 10 3 kg/m 3 , 
and a layer of Au of thickness A^ u — 203 nm. The plate of density p$i = 2.33 x 10 3 kg/m 3 
was coated first with the same thickness of Cr and then with a layer of Cu of thickness 
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Acu = 200 nm. Considering that the conditions z, A <C R, D are satisfied, the hypothetical 
force is given by j^, 1^ 



F*»{z) = -A-K 2 Ga\ 3 e-*/ x R [p Au - (p Au - p Cr ) e - A ^ x - ( PCr - p.) e -^ +A ^ x ] (31) 

x [pen ~ (Pen - Per) e~ A ^ x - (per - P Si) e^+^A] . 

In our experiment, the strongest constraints on the Yukawa hypothetical interactions are 
obtained from the dynamic measurement of the parallel plate pressure (i.e., Casimir force 
gradient) rather than the static measurement of the Casimir force. For this case, the Yukawa 
pressure can be found from Eq. (|31|) by the using Eq. which follows from the proximity 
force theorem, Eq. (jSj), 

phy P{z) = _2nGa\ 2 e- z / x [p Au - (p Au - p Cr ) e~ A ^ /x - (p Cr ~ Ps) e~^ +A ^ x ] (32) 

X [pen ~ (PCu - PCr) e' A ^' X - (p C r ~ P Si) e'^ +A ^' X ] . 

Note that the Newtonian gravitational pressure is also below the sensitivity of the present 
experiment and can be neglected. 

As shown in Ref. j^], surface roughness can significantly influence the magnitude of a 
hypothetical force in the nanometer range. To compute the hypothetical pressure taking 
account of roughness, one can use exactly the same method that was applied in Sec. IV in 
the case of the Casimir pressure. The result is 

n 

Pr P (z) = E v i v 3 phyP (* + 2H ° - h i- hj), (33) 

where our notation is that of Sec. IV, and P hyp is given by Eq. ()32j). 

With these results, we can now obtain constraints on the hypothetical Yukawa pressure 
from the agreement of our measurements of the Casimir pressure with theory. According to 
the results of Sec. V, the optimal separation region for obtaining constraints is z > 310 nm. 
Here the r.m.s. deviation between theory and experiment, o"^ 00 = 0.34 mPa, is somewhat 
smaller than would be the case if the complete measurement range were used, which includes 
the shortest separations 260 nm< z <310nm. (At the shortest separations the experimental 
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relative error is less than the theoretical error, but it increases with separation so it even- 
tually exceeds the theoretical relative error.) Within the interval z > 310 nm the strongest 
constraints are obtained from the shortest separations. We choose z = 320 nm and obtain 
constraints from the inequality 



In Fig. El constraints on a following from Eq. (|34|) are plotted for different values of 
the interaction range A (curve 1). In the same figure constraints from previous experiments 
are also shown. They were obtained from old measurements of the Casimir force between 
dielectrics js^j (curve 2), from Casimir force measurements by means of a torsion pendulum 



region in the (a, A) plane above the curve is excluded, and below the curve is allowed 
by the experimental results. As can be seen from Fig. El the present experiment leads 
to the strongest constraints in a wide interaction range, 56nm< A <330nm. The largest 
improvement, by a factor of 11, is achieved at A ~ 150 nm. We note that the constraints 
obtained here almost completely fill in the gap between the modern constraints obtained 
by AFM measurements, and those obtained using a torsion pendulum. Within this gap 
the best previous constraints were obtained from old measurements of the Casimir force 
between dielectrics which are not as precise or reliable as those obtained here from the 
Casimir pressure measurements between metals using an MTO. 

Turning to the power- law- type hypothetical interactions given by Eq. (|29j) . the present 
experiment does not lead to improved constraints. This is explained by the fact that the 
metallic coatings used were too thin (and hence, too light) to give a significant contribution to 
hypothetical interactions with a longer interaction range. In fact, the thicker (bulk) matter 
contributes more significantly in this case, even though its density is much lower than for the 
metal coatings. To obtain stronger constraints on the constants characterizing new power- 
law interactions, thicker metal coatings and larger interacting bodies are preferable. We 
anticipate that future measurements of the Casimir force will employ such samples. 




(34) 




AFM 




(curve 4). In all cases the 
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VII. CONCLUSIONS AND DISCUSSION 



Our primary objective in the present paper has been to set new limits on extra- 
dimensional models and other physics beyond the Standard Model using Casimir force 
measurements between a sphere and a plate separated by ~ 0.2-1.2 /im. These experi- 
mental results along with a detailed theoretical analysis lead to new constraints on Yukawa 
modifications of Newtonian gravity at short distances, and these are presented in Fig. 14. 

Although the constraints on the Yukawa parameter a we obtain are ~ 10 13 , this does 
not imply that our experiment has to be improved by 13 orders of magnitude to detect 
Newtonian gravity. Rather, the exponential factor e~ r l x in Eq. (J26|) strongly suppresses 
contributions from those parts of the sphere and plate separated by > A. This "finite size 
effect" , which is not relevant for gravity, has a consequence that a Yukawa force between 
the sphere and plate with a = 1 is much weaker than gravity. The actual gravitational force 
in our experiment is only ~ 5 orders of magnitude from being detected (using the actual 
materials employed), a gap that may be closed within the forseeable future based on the 
rapid progress various groups have made in the last few years. 

To carry out our objective, an MTO was used to obtain the first precise measurements of 
the Casimir force between dissimilar metals Cu and Au. In the static regime, the Casimir 
force between a sphere and a plate was measured with an absolute error 0.3 pN at a 95% 
confidence level. This translates into an experimental relative error of the Casimir force 



times smaller 
21. 
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measurements at the separation 188 nm of approximately 0.27%, i.e., severa 
than the most precise previous experiments carried out by means of an AFM 
To take advantage of the high quality factor of the MTO, the Casimir force derivative 
between the sphere and plate was then measured dynamically. This derivative was shown to 
be effectively equivalent to the pressure between two parallel plates composed of the same 
materials. This pressure was determined with a mean absolute error of 0.6 mPa at the same 
confidence level, which leads to an experimental relative error of approximately 0.26% at 
the separation 260 nm. (This compares with a relative error of 15% quoted by the authors 



of Ref. 



26| who directly measured the Casimir pressure between parallel plates.) 



As noted above, the precise calculation of the Casimir force and pressure between real 
metals that is needed to extract the effects of new physics calls for the careful accounting of 
different corrections (e.g., due to surface roughness, finite conductivity of the boundary met- 
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als, nonzero temperature, and finite extent of the plate). The corrections due to roughness 
were computed on the basis of the AFM images of the surfaces of the test bodies. Finite 
conductivity corrections were calculated using tabulated data for the complex index of re- 
fraction. Other corrections were also estimated (including traditional thermal corrections) 
and found to be negligible. The estimated theoretical uncertainty is at the level of 1% of the 
calculated force, i.e., greater than the experimental uncertainty at the shortest separations 
used. Within the limits of all errors, theory is in good agreement with experiment. 

The level of agreement between theory and experiment was used to draw important 
conclusions concerning the influence of thermal effects on the Casimir force predicted by 
quantum field theory at nonzero temperature. Our experimental results lead to a resolution 
of the controversy over whether the thermal effects on the Casimir force between real metals 
is close to that predicted by the Matsubara formulation of quantum field theory for ideal 



metals, or is significantly different as claimed in Refs. 



42, 43, 441. We have shown that 



the experimental results contradict the large, linear in temperature, thermal corrections 



predicted in Refs. |42j, |4 



44|. Although the sensitivity of the current experiment is not yet 



sufficient to detect the small thermal corrections to the Casimir force and pressure predicted 
for real metals in Refs. 0, 33^ they are compatible with our experiment. 

The good agreement between theory and experiment was then used to set stronger con- 
straints on hypothetical Yukawa interactions predicted by extra-dimensional physics and 
extensions to the Standard Model. Existing limits were strengthened by a factor of up to 
11 within a wide interaction range, from 56 nm to 330 nm. This interaction range covers the 
gap between the modern results obtained from Casimir force measurements using a torsion 
pendulum and an AFM. The previously known constraints within this gap were based on 
old, and less reliable, Casimir force measurements between dielectrics. 

Our experimental arrangement suggests the need for additional work to further improve 
the agreement between theory and experiment. We plan to use smoother and thicker metal 
coatings on the surfaces of the test bodies, yielding a better characterization of the rough- 
ness. This will permit us to carry out measurements at shorter separations, and to sig- 
nificantly strengthen the constraints on the predictions of extra-dimensional physics in a 
wider interaction range. The altimate goal of this program to use the iso-electronic effect 



18, 



39 



78, 



79 



80] to suppress the Casimir force so as to improve our sensitivity to new 



forces beyond the Standard Model. 
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Figures 



FIG. 1: Schematic of the experimental setup showing its main components, see text. Inset: Res- 
onance curve for the MTO. Also shown is a Lorentzian fit with Q ~ 8000. See text for further 
details. 

FIG. 2: Schematic of the bridge circuit used to measure the capacitance. The DC voltages V\ and 
V2 are used to correct for small deviations when no interactions are present. They also linearize 
the response of the circuit. Details of the charge amplifier and the part of the circuit to balance 
the bridge are omitted for clarity. Inset: Scanning electron microscopy image of the MTO. 

FIG. 3: Dependence of the angular deviation 9 as a function of the applied voltage to the sphere. 
Data obtained at two different separations z between the metallic layers are shown. 

FIG. 4: Electrostatic force F e as a function of separation z for AV = Va u — = 0.27 V and 
AV = 0.35 V. The fits of the data using Eq. ((2j) are also shown by the solid lines. See text for 
further discussion. 

FIG. 5: Angular displacement of the MTO as a function of linear displacement when the sphere is 
moved parallel to the MTO's axis (top axis, A) and when it is moved perpendicular to the MTO's 
axis (bottom axis, v)- The data were acquired at z = 3 /im with AV = 0.27 V. For comparison, 
the expected values for an infinite Cu layer using Eq. © are shown as solid lines. 

FIG. 6: lxl /mi 2 atomic force microscopy image of the Cu layer. The topography of the Au layer 
on the sphere is similar. The gray scale to the right of the image gives the height of peaks above 
the bottom of roughness, with lighter tones corresponding to larger heights. The set of such images 
was used to analyze the effect of the surface roughness on the Casimir force as described in the 
text. 
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FIG. 7: Relief heights hi versus the fraction of the total area with heights h < hi+\. 



FIG. 8: Absolute value of the measured Casimir force as a function of separation obtained using 
the static mode. The value of the separation between the two metals is determined as discussed in 
the text. 



FIG. 9: Absolute value of the parallel plate Casimir pressure as a function of separation obtained 
from the dynamic measurement. The value of the separation between the two metals is determined 
as discussed in the text. 



FIG. 10: Difference of the theoretical and experimental Casimir forces between the sphere and 
plate versus separation obtained from the static measurement. 



FIG. 11: Difference of the theoretical and experimental parallel plate Casimir pressures versus 
separation obtained from the dynamic measurement. 
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FIG. 12: Difference of the theoretical parallel plate Casimir pressure as predicted by Refs. 
(which incorporates an alternative thermal correction) and experiment versus separation. 



FIG. 13: Difference of the theoretical parallel plate Casimir pressures, as predicted by Ref. 
(which incorporates another alternative thermal correction) and experiment versus separation 



44] 



FIG. 14: Constraints on the Yukawa interaction constant a versus interaction range A. Curve 
1 is obtained in this paper, curve 2 follows from old measurements of the Casimir force between 
dielectrics 35|. Curves 3 and 4 are obtained from the Casimir force measurements between metals 



by use of the torsion pendulum h]j 3 an d by means of an atomic force microscope 



3Q, 



respectively. The region in the (a, A) plane above each curve is excluded, and below each curve is 
allowed. 
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